JOURNAL OF VIROLOGY, May 2005, p. 6180-6193 
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.10.6180-6193.2005 


Vol. 79, No. 10 


Copyright © 2005, American Society for Microbiology. All Rights Reserved. 


Comparative Host Gene Transcription by Microarray Analysis Early 
after Infection of the Huh7 Cell Line by Severe Acute Respiratory 
Syndrome Coronavirus and Human Coronavirus 229E 


Bone S. F. Tang,’ Kwok-hung Chan,' Vincent C. C. Cheng,' Patrick C. Y. Woo,’ Susanna K. P. Lau,’ 
Clarence C. K. Lam,” Tsun-leung Chan,’ Alan K. L. Wu,’ Ivan F. N. Hung,' Suet-yi Leung,” 
and Kwok-yung Yuen'* 


Department of Microbiology, Centre of Infection and Immunology, State Key Laboratory of Emerging Infectious Diseases, 
The University of Hong Kong,‘ and Division of Haematology” and Division of Anatomical Pathology,? Department of 


Pathology, The University of Hong Kong, Hong Kong 


Received 21 October 2004/Accepted 11 January 2005 


The pathogenesis of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) at the cellular 
level is unclear. No human cell line was previously known to be susceptible to both SARS-CoV and other human 
coronaviruses. Huh7 cells were found to be susceptible to both SARS-CoV, associated with SARS, and human 
coronavirus 229E (HCoV-229E), usually associated with the common cold. Highly lytic and productive rates of 
infections within 48 h of inoculation were reproducible with both viruses. The early transcriptional profiles of 
host cell response to both types of infection at 2 and 4 h postinoculation were determined by using the 
Affymetrix HG-U133A microarray (about 22,000 genes). Much more perturbation of cellular gene transcription 
was observed after infection by SARS-CoV than after infection by HCoV-229E. Besides the upregulation of 
genes associated with apoptosis, which was exactly opposite to the previously reported effect of SARS-CoV in 
a colonic carcinoma cell line, genes related to inflammation, stress response, and procoagulation were also 
upregulated. These findings were confirmed by semiquantitative reverse transcription-PCR, reverse transcrip- 
tion-quantitative PCR for mRNA of genes, and immunoassays for some encoded proteins. These transcrip- 
tomal changes are compatible with the histological changes of pulmonary vasculitis and microvascular 


thrombosis in addition to the diffuse alveolar damage involving the pneumocytes. 


Severe acute respiratory syndrome (SARS)-associated coro- 
navirus (SARS-CoV) is the etiological agent of SARS (26, 39). 
The disease is associated with significant mortality and mor- 
bidity (38). Such aggressive clinical behavior is very different 
from that of other known human coronaviruses, such as the 
group 1 coronaviruses 229E and NL63 or the group 2 corona- 
virus OC43. Although these viruses are generally associated 
with mild upper respiratory tract infection such as the common 
cold (47-49), these human coronaviruses can also cause pneu- 
monia when the very young, elderly, or immunosuppressed 
hosts are affected (14, 15, 40). Animal models of SARS were 
established by infecting cynomolgus monkeys, ferrets, domestic 
cats, or mice (16, 27, 31, 33, 55), but there are questions about 
their reproducibility of the pathology mimicking human dis- 
ease. Although the cellular receptors for attachment of SARS- 
CoV were found to be ACE2 (30) and, recently, L-SIGN (23), 
pathogenesis at the cellular level is largely unknown. 

SARS-CoV is an enveloped, positive-sense, single-stranded 
RNA virus which can grow in embryonal monkey cell lines, 
including Vero E6 and fetal rhesus monkey kidney (FRhk-4) 
cells (26, 39). It can be subcultured onto other Vero cells and 
colonic carcinoma cell lines such as Caco-2 or LoVo (3, 10). 
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Unlike other human coronaviruses, SARS-CoV proliferates 
rapidly and causes obvious cytopathic effects in Vero E6 cells 
within 48 h of inoculation (35). There are no other human cell 
lines known to be susceptible to infection by both SARS-CoV 
and other human coronaviruses. Recently there have been 
reports of a human hepatoma cell line, Huh7, which can be 
infected by pseudotyped lentiviral particles carrying the spike 
protein of SARS-CoV and wild-type replicative SARS-CoV (6, 
20, 45). We report in this study the susceptibility of the cell line 
Huh7 to infection by both SARS-CoV and human coronavirus 
229E (HCoV-229E). A comparative gene transcriptional pro- 
file at an early stage of infection of Huh7 cells by these two 
viruses was performed to elucidate their difference and its 
importance in the pathogenesis of disease. 


MATERIALS AND METHODS 


Cell lines and virus. The human hepatoma cell line Huh7 (courtesy of David 
Ho, Aaron Diamond AIDS Research Center) was used throughout this study. 
The cells were incubated at 37°C in minimal essential medium supplemented 
with 10% fetal calf serum and 100 IU/ml of penicillin and 100 g/ml of strep- 
tomycin. Our prototype virus (SARS-CoV HKU-39849) was isolated from the 
lung tissue biopsy of the brother-in-law of the index patient who traveled from 
Guangzhou and started a superspreading event in the Hong Kong Special Ad- 
ministrative Region leading to the pandemic (37). The HCoV-229E strain 
(ATCC VR-740) was used in this study. The SARS-CoV and HCoV-229E strains 
used in our experiments had undergone three passages in the FRhk-4 cell and 
MRC-5 cell lines, respectively, and were stored at —70°C. Viral titers were 
determined as median tissue culture infective dose (TCIDs 9) per ml in confluent 
Huh7 cells in 96-well microtiter plates, which standardized the viral inoculum 
and measured the relative susceptibility of the Huh7 cell line to these two viruses. 
The relative susceptibilities of the Vero 1008, Vero 76, Vero, and Huh7 cell lines 
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to SARS-CoV and HCoV-229E were also tested by TCIDs59. One hundred 
TCIDs9 was confirmed by plaque assays to be equivalent to 85 PFU. All work 
with infectious virus was performed inside a type II biosafety cabinet in a 
biosafety containment level III facility, and the personnel wore powered air- 
purifying respirators (HEPA Airmate; 3 M, Saint Paul, Minn.). 

Monitoring of virus-induced cytopathic effect, antigen detection, and semi- 
quantitative and quantitative RT-PCR. Huh7 cells and culture supernatants 
infected with either SARS-CoV or HCoV-229E at a multiplicity of infection of 
100 TCIDsp per cell were collected at 2, 4, 12, and 24 h postinfection. A washing 
step was performed 1 h postinoculation. The percentages of cells developing 
cytopathic effects were counted by inverted light microscopy at 24 and 48 h. The 
rate of viral replication was measured by reverse transcription-quantitative PCR 
(RT-qPCR) on the culture filtrate. The amount of coronavirus antigen expres- 
sion in infected cells was measured by indirect immunofluorescence with conva- 
lescent-phase sera of patients with SARS-CoV or HCoV-229E infection, as 
reported previously (2, 8). Briefly, harvested cells were prepared and fixed in 
ice-cold acetone for 10 min. Convalescent-phase serum at a dilution of 1 in 100 
was used to react with the infected cells harvested at various time points. After 
30 min incubation, the cells were washed twice in phosphate-buffered saline for 
5 min each, and then goat anti-human fluorescein isothiocyanate conjugate 
(INOVA Diagnostics, Inc., San Diego, CA) was added and the cells were further 
incubated for 30 min at 37°C. The cells were washed again as described above, 
and the percentage of positive cells was manually estimated under UV micros- 
copy. 

Reverse transcription-PCR (RT-PCR) for SARS-CoV and HCoV-229E was 
done directly on culture filtrate according to our previous protocol (2). Briefly, 
total RNA extracted from culture filtrate with the QiAamp virus RNA minikit 
(Qiagen) as instructed by the manufacturer was reverse transcribed with random 
hexamers. cDNA was amplified with SARS-CoV primers (forward, 5’-TACAC 
ACCTCAGCGTTIG-3’; reverse, 5’-CACGAACGTGACGAAT-3’) and HCoV- 
229E primers (forward, 5'-GGTACTCCTAAGCCTTCTTCG-3’; reverse, 5'-G 
ACTATCAAACAGCATAGCAGC-3’). Using real-time RT-qPCR assays, 
cDNA was amplified in SYBR Green I fluorescence reactions (Roche, Mann- 
heim, Germany). For the RT-qPCR of SARS-CoV, a 20-yl reaction mixtures 
containing 2 pl cDNA, 3.5 mmol/liter magnesium chloride, and 0.25 pmol/liter 
of the same forward and reverse primers as in the reaction mixtures were thermal 
cycled with a Light Cycler (Roche) (95°C for 10 min, followed by 50 cycles of 
95°C for 10 s, 57°C for 5 s, and 72°C for 9 s with ramp rates of 20°C/s). For the 
RT-qPCR of HCoV-229E, the conditions were similar to those described above 
except that 45 cycles of 95°C for 10 s, 65°C for 3 s, and 72°C for 12 s were used. 
Plasmids with the target sequences were used to generate the standard curve. At 
the end of the assay, the PCR products of SARS-CoV and HCoV-229E (182 bp 
and 295 bp, respectively) were subjected to a melting curve analysis (65 to 95°C, 
0.1°C/s) to determine the specificity of the assay. All assays were performed in 
replicates. 

Microarray analysis. Human genome-wide gene expression was examined 
with the GeneChip system HG-U133A microarray (Affymetrix Inc., Santa Clara, 
CA), which is composed of more than 22,000 oligonucleotide probe sets inter- 
rogating approximately 18,400 unique transcripts, including 14,500 well-charac- 
terized human genes. Quality control, GeneChip hybridization, and data acqui- 
sition and analysis were performed at the Genome Research Centre, The 
University of Hong Kong, according to the standard protocols available from 
Affymetrix. In brief, total RNAs of the infected or uninfected cell lines at 
different time points were extracted using the RNeasy minikit (Qiagen, Valencia, 
CA). Double-stranded cDNA was synthesized from 10 wg of total RNA with the 
GeneChipT7-Oligo (dT) Promoter Primer Kit (Affymetrix, Inc.) and the Super- 
Script Choice System (Invitrogen). Biotin-labeled cRNA was then synthesized by 
in vitro transcription using the BioArray High Yield RNA Transcript Labeling 
Kit (Affymetrix, Inc). After fragmentation, 15 jg of labeled CRNA was hybrid- 
ized to the oligonucleotide microarray. The chips were washed and stained using 
the GeneChip Fluidics Station 400 (Affymetrix) and then scanned with the 
GeneChip Scanner 3000 (Affymetrix). Data analysis was performed using the 
Microarray Suite Expression Analysis software (version 5.1; Affymetrix). For 
comparison across different arrays, the data for each array were normalized by a 
global scaling strategy, using a scaling target intensity of 500. By using the 
Affymetrix-defined comparison mathematical algorithms, a fold change in ex- 
pression between each of the infected samples in comparison to the uninfected 
mock control was calculated, log, transformed, and further classified as not 
changed, increased (signal log ratio change P value of <0.005), decreased (signal 
log ratio change P value of >0.995), or marginally increased or decreased. To 
classify a gene as significantly upregulated or downregulated after infection at a 
specific time point, two additional criteria were used: (i) the fold change must be 
greater than or equal to 2 (signal log ratio of =1 if upregulated or =1 if 
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downregulated) to be classified as increased or decreased, and (ii) genes that 
were classified as upregulated must be flagged as present in the infected samples, 
while genes that were classified as downregulated must be flagged as present in 
the uninfected control sample. All gene chip procedures were performed in 
replicates. 

Gene expression analysis by semiquantitative PCR, RT-qPCR, and immuno- 
assays. Genes with significant transcriptional changes known to be associated 
with biological significance were selected for further analysis by semiquantitative 
PCR, RT-qPCR, and immunoassays. RT-qPCR was performed according to our 
previous protocol (56). The extracted RNA was pretreated with DNase. Primers 
which specifically amplified nine genes related to apoptosis, inflammation, and 
coagulation were designed (Table 1). First-strand cDNA was synthesized from 
the total RNA by reverse transcription with random hexamers. Semiqualitative 
comparison was performed using simple gel electrophoresis and ethidium bro- 
mide staining (10). Quantitative PCR was performed using the SYBR Green I 
fluorescence reactions in a Light Cycler as described above. Detailed PCR 
conditions are available upon request. Serial dilutions of a reference cDNA 
derived from the SARS-CoV-infected sample were used to generate the standard 
curve. Melting curve analysis was performed for each primer pair at the end of 
the reaction to confirm the specificity of the assay. The housekeeping porpho- 
bilinogen deaminase gene was used for standardization of the initial RNA con- 
tent of a sample. Experiments were performed in duplicate, and the result for an 
individual sample was expressed as the mean expression level of a specific 
gene/porphobilinogen deaminase gene relative to the reference cDNA. The 
relative expression between each infected sample and the uninfected control was 
then calculated and expressed as fold change. 

Three sets of immunoassays (human interleukin-8 [IL-8] from BD Biosciences 
and PAI1 [serpine 1] and TFPI2 from Diagnostica Stago) were performed 
according to our previous protocols (5, 28) and the manufacturers’ instructions. 
Briefly, for IL-8 assay, 100 jl culture supernatant or standard was added to IL-8 
monoclonal antibody-coated wells and incubated for 2 h at room temperature. 
After washing five times, 100 pl a biotinylated anti-human IL-8 monoclonal 
antibody was added to each well. After five washings, 100 yl of tetramethyl 
benzidine substrate was added, followed by incubation for 30 min at room 
temperature. The reaction was stopped by adding 50 wl stop solution. The 
absorbance was read at 450 nm within 30 min. The experimental conditions and 
procedure were similar for the immunoassays for PAI1 and TFP1. All experi- 
ments were done in replicates to ensure accuracy and reproducibility. 

Statistical analysis. The fold change of the target gene expressions and the 
difference in the concentrations of expressed proteins between SARS-CoV and 
HCoV-229E at different postinoculation time points were compared by Student’s 
t test. A P value of <0.05 was considered significant. A statistical package (SPSS 
10.0; SPSS Hong Kong, Hong Kong) was used for all analyses. 


RESULTS 


Susceptibility of the Huh7 cell line to SARS-CoV and HCoV- 
229K. With a multiplicity of infection of 100, a cytopathic effect 
was visible in Huh7 cells at 24 and progressed to about 50% 
cell death at 48 h with both viruses. SARS-CoV but not HCoV- 
229E grew well in the Vero cell line. The viruses produced 
comparable TCIDs9 of around 107 per ml in the culture su- 
pernatant of Huh7 cells at 48 h (Table 2). In terms of viral 
load, a 1-log-unit increase of viral genome copy was noted at 
12 h postinoculation with both viruses, which was followed by 
a peak at 24 h (Table 3). Antigen expression could be observed 
by indirect immunofluorescence in over 50% of the cells at 
24 h postinoculation for both viruses. 

Effects on gene expression of host cells determined by mi- 
croarray. Two hundred twenty-four genes were significantly 
altered within 4 h postinoculation in the transcriptomal expres- 
sion analysis. Only 21 genes were perturbed by HCoV-229E 
per se, whereas 164 genes were altered by SARS-CoV infec- 
tion only, and the remaining 39 genes were altered by both 
coronaviruses. Out of the 164 genes with altered expression in 
SARS-CoV, 38 were upregulated and only one was downregu- 
lated at both 2 and 4 h postinoculation. Forty-three were up- 
regulated and 16 were downregulated at 2 h postinoculation. 
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TABLE 1. Primers used for checking the expression of nine selected genes and a house keeping gene by semiquantitative RT- PCR 


and RT-qPCR 
Gene product Primer Sequence (5’ to 3’) ee _ pene et 

Tissue factor pathway inhibitor 2 TFPI2-F CACCTATACTGGCTGTGGAGG 21 236 
TFPI2-R CCTCATGCTGTCATATTATTCTTC 24 

Plasminogen activator inhibitor 1 (serpine 1) PAII-F CCGAGGAGATCATCATGGAC 20 259 
PAII-R GCCAAGGTCTTGGAGACAGA 20 

Thrombospondin 1 THBS1-F ATCATGGCTGACTCAGGACC 20 229 
THBS1-R GGAAGCCAAGGAGAAGTGAT 20 

Interleukin-8 IL8-F TGGAAGAGAGCTCTGTCTGGA 21 215 
IL8-R CCAGGAATCTTGTATTGCATCT 22 

Nuclear factor of kappa light polypeptide gene NFKB2-F CCTGGCAGGTCTACTGGAGG 20 281 
enhancer in B-cells 2 NFKB2-R AAATAGGTGGGGACGCTGT 19 

Jun B proto-oncogene JUNB-F CCAGCTCAAACAGAAGGTCATG 22 209 
JUNB-R GGAGTCCAGTGTGGTTTGCG 20 

Pleckstrin homology-like domain, family A, PHLDAI-F CTCTCATCCTCACTCGCACC 20 214 
member 1 PHLDA1-R TTGATCCAAGTGAAGACAATAGAA 24 

Caspase recruitment domain family, member 10 CARD10-F CACGTGGAGGTGACTGAGAA 20 274 
CARD10-R CTCAGGCCTCACTGCTGCT 19 

BCL2-associated X protein BAX-F GATGATTGCCGCCGTGGACA 20 175 
BAX-R GATGGTTCTGATCAGTTCCGG 21 

Porphobilinogen deaminase PBGD-F AAGATGAGAGTGATTCGCGTG 21 226 
PBGD-R GGTCCACTTCATTCTTCTCCAG 22 


Forty-nine were upregulated and 17 were downregulated at 4 h 
postinoculation. In contrast, for HCoV-229E infection, only 
one gene was upregulated and no gene was downregulated at 
both 2 and 4 h postinoculation. No gene was upregulated and 
two were downregulated at 2 h postinoculation. Fourteen 
genes were upregulated and four were downregulated at 4 h 
postinoculation. When multiple transcripts of the same gene 
were eliminated and analyzed, genes related to apoptosis (n = 
23), inflammatory or immune response (n = 34), and coagu- 
lation (n = 5) were identified, in addition to the expected 
genes related to the stress response and metabolism and other 
unknown genes (Table 4). Of the 23 apoptotic genes affected, 
13 of them are proapoptotic and 11 of them are upregulated in 
SARS-CoV infection, compared with only 3 in HCoV-229E 
infection (Table 5). For inflammation and immune response, 
32 genes are upregulated in SARS-CoV, compared with only 3 
in HCoV-229E. These include the genes for NFKB1A; 


TABLE 2. Viral load of SARS-CoV and HCoV-229E in different 
cell lines 48 h postinoculation 


Logig TCIDs59/ml 


Cell line 
SARS Co-V HCoV-229E 
Vero 1008 6.5 ND* 
Vero 76 8.33 ND 
Vero 6.5 ND 
Huh7 7.0 6.5 


“ND, not detectable. 


NFKB2; IL-8; transforming growth factor B2; chemokines 
CXCLI, -2, -3, -5, -6, and -10; ICAM1; and tumor necrosis 
factor alpha-induced proteins. Quite unexpectedly, genes of 
the procoagulation pathway were also affected by SARS-CoV 
infection, with upregulation of PLSCR1 (phospholipid scram- 
blase 1), EGR1 (early growth response 1 gene), PAT1/SER- 
PINE1 (plasminogen activator inhibitor 1), and THBS1 
(thrombospondin 1). In terms of stress response, seven genes 
were upregulated in SARS-CoV infection, compared with only 
one in HCoV-229E infection. Overall there were far more 
changes in expression of genes related to the cell cycle, tran- 
scription, and metabolism and those with miscellaneous and 
unknown functions in SARS-CoV infection. When the Path- 
way Assist software (Ariadne Genomics Inc.) was used for 
linking the altered genes in cellular pathways for SARS-CoV, 
there were clear clusterings of altered genes related to apo- 
ptosis, inflammation, and coagulation (Fig. 1). 


TABLE 3. Number of viral genome copies of SARS-CoV and 
HCoV-229E at different postinoculation time points in the Huh7 


cell line 
Copies/ml 
h postinoculation 

SARS-CoV HCoV-229E 

2 2.2 X 107 3.8 x 107 

4 1.5 x 107 3.3 x 10” 

12 3.6 x 10° 7.8 x 108 

24 1.2 x 10° 3.4 x 10'° 
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TABLE 4. Numbers of genes with significant change of expression 
in different categories and subcategories 


No. of genes 


Category and 


subcategory Total SARS- HCoV- Both SARS-CoV 
oe Cov 229 and HCoV-229E 

Apoptosis 23 

Antiapoptosis 3 3 0 0 

Proapoptosis 12 7 2 3 

DNA repair 5 2 2 1 

Other 3 2 1 0 
Inflammation 34 

Antiinflammatory 2 2 0 0 

Proinflammatory 24 21 0 3 

Immunomodulatory 2 1 1 0 

Other 6 5 1 0 
Coagulation 5 

Anticoagulants 1 0 0 1 

Procoagulants 4 4 0 0 
Metabolism 35 

Amino acid 4 2 1 il 

Carbohydrate 7 6 0 1 

Lipid a} =) 0 0 

Transporter 10 6 1 3 

Other 9 7 2 0 
Signal Transduction 18 15 1 2 
Transcription 37 25 3 9 
Cell Cycle 17 12 3 2 
Cell growth and ¢ 6 0 1 

differentiation 
Stress response 7 6 0 1 
Miscellaneous 16 11 1 4 
Unknown 25 16 2 7 
Total 224 164 21 39 


Confirmation of cellular gene and protein expression by 
semiquantitative PCR, RT-qPCR, and immunoassays. Since 
previous studies have documented the alterations of the apo- 
ptotic and inflammatory pathways in relation to SARS-CoV 
infection, only key genes related to inflammation and apopto- 
sis were chosen to confirm the results of microarray analysis by 
semiquantitative RT-PCR and RT-qPCR. For the coagulation 
pathway, both procoagulation and anticoagulation genes were 
chosen. A similar trend of upregulation of gene expression was 
found by RT-qPCR, with SARS-CoV showing a 1.4- to 10.8- 
fold increase compared with HCoV-229E for coagulation 
(TFPI2, PAI1, and THBS1), inflammation (IL-8 and NFKB2), 
transcription (JUNB), and apoptotic (PHLDA1, CARD10, 
and BAX) pathways (Fig. 2 and 3). IL-8 was increased in both 
SARS-CoV and HCoV-229E infection, but the increase was 
much higher in the case of SARS-CoV. The enzyme immuno- 
assay showed that SARS-CoV induced higher concentrations 
of PAI1 and IL-8 than HCoV-229E at 2, 4, 12, and 24 h 
postinoculation. Both SARS-CoV and HCoV-229E induced 
similar TFPI2 expression at 4, 12, and 24 h postinoculation, but 
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at 2 h postinoculation SARS-CoV induced a lower concentra- 
tion of this predominately anticoagulation protein (Fig. 4). 


DISCUSSION 


Of the four coronaviruses known to infect human, HCoV- 
229E, HCoV-OC43, and NL63 are generally associated with 
mild upper respiratory tract infection such as the common cold 
in immunocompetent hosts (47-49). In contrast, SARS-CoV 
causes respiratory failure in over 60% of affected persons, with 
a mortality rate of 15% (37, 38). Besides pneumonia, SARS is 
also clinically manifested as watery diarrhea without enteroco- 
litis (48.6%) (8), hepatitis without liver failure (49.4%) (22), 
lymphopenia (75%) (37), impaired coagulation (63%) (51), 
and occasionally pulmonary vasculitis and thrombosis in the 
lungs of those who died (36, 13). Much is already known about 
many aspects of SARS, including the virology (18, 26, 39, 53), 
genomics (32), diagnostics (5, 28, 43, 52, 54), clinical features 
and progression in relation to viral load (22, 37), treatment (7, 
9, 46), infection control (44), and immunization (55). However 
little is known about the pathogenesis at the cellular level, 
despite the identification of ACE2 and L-SIGN as the recep- 
tors for binding of SARS-CoV to host cells (30, 23). 

One important limitation to studies at the cellular level is the 
absence of a relevant cell line which can be lytically infected by 
SARS-CoV. Initially, only embryonal monkey kidney cell lines 
such as Vero or FRhk-4 and their derivatives could be readily 
infected. Subsequent investigation showed that human colon 
adenocarcinoma cell lines Caco-2, CL-14, and LoVo could also 
be infected (3, 10). However, inflammation or cell death is not 
manifested on endoscopic examination of the intestines of 
SARS patients (29). Examination of intestinal tissue biopsy by 
electron microscopy revealed abundant intracellular viral par- 
ticles, and there was negligible inflammatory cells or cellular 
apoptosis on light microscopy (29). In fact, a high TCID.p in 
the culture supernatant of persistently infected LoVo cells can 
be achieved without any cytopathic effects (3). Recently a gene 
expression analysis showed that SARS-CoV upregulated anti- 
apoptotic genes and several CXC chemokines while downregu- 
lating proapoptotic genes, IL-18, and macrophage migration 
inhibitory factor (10) at 24 h postinoculation of Caco-2 cell 
lines. These patterns of gene expression appeared to explain 
the severe watery diarrhea without clinical or endoscopic signs 
of necrosis or inflammation in SARS patients (29). However, 
no data on early gene expression before the onset of cytolysis 
were reported. Similarly, no comparative study of the host cell 
transcriptional profile and those for other, less virulent coro- 
naviruses was ever reported. 

No pneumocyte cell line has yet been found to support lytic 
or nonlytic infection by SARS-CoV. However, the human hep- 
atoma cell line Huh7 was found to be susceptible to SARS- 
CoV. This is not completely unexpected, because the mouse 
hepatitis virus, a group 2 coronavirus, is known to infect Huh7 
cells (25). In this study HCoV-229E was found to produce 
good lytic infection within 48 h postinoculation. A high multi- 
plicity of infection of 100 TCID.,s per cell was used to ensure 
the reproducibility of gene expression, as previously reported 
(41). Since the expression of a high number of genes was 
expected to change significantly when virus-induced cytopa- 
thology was impending for a rapidly lytic viral infection, we 
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TABLE 5. Signal log ratios (log, scale) of the genes that are significantly altered after SARS-CoV or HCoV-229E infection 


Signal log ratio (log, scale)” 


Ene after infection with: 
Category ane F Description 
subcategory Affymetrix P HCoV-229E SARS-CoV 
identification Designation 
no. 2h 4h 2h 4h 
Apoptosis 
Antiapoptosis 200797_s_at MCL1 Myeloid cell leukemia sequence 1 (BCL2 related) 0.2 0 1 0.3 
210538_s_at BIRC3 Baculoviral IAP repeat containing 3 0.2 0.6 2.1 14 
201631_s at ITER3 Immediate-early response 3 0.3 —0.2 L7 1:3, 
Apoptosis 204285_s_at PMAIP1? Phorbol-12-myristate-13-acetate-induced protein 1 0.5  —0.2 1.6 0.6 
204908 _s at BCL3? B-cell CLL/lymphoma 3 0.1 0.1 1.9 1. 
207181_s_ at CASP7 Caspase 7, apoptosis-related cysteine protease —0.2 0.4 0.1 12 
208478_s_ at BAX BCL2-associated X protein 0.5 0.6 0.4 11 
208536_s_at BCL2L11 BCL2-like 11 (apoptosis facilitator) e135 0.3 0.1 12 
209305_s at GADD45B? Growth arrest and DNA-damage-inducible, beta 0.1 —0.2 14 0.8 
210025_s at CARD10? Caspase recruitment domain family, member 10 1 —0.4 19 2.7 
211085_s_ at STK4 Serine/threonine kinase 4 0.4 1:3 0.2 0.6 
211695_x_at MUC1 Mucin 1, transmembrane O02. e135: 0.1 0.8 
218368_s_ at TNFRSF12A Tumor necrosis factor receptor superfamily, 0.5 0 i) 0.6 
member 12A 
218880_at FOSL2 FOS-like antigen 2 18 —0.4 3.3 15 
217997 at? PHLDA1 Pleckstrin homology-like domain, family A, iS 0.1 2.9 1.6 
member 1 
DNA repair 205022_s_at CHES1 Checkpoint suppressor 1 0 —0.2 0.2 =1 
205072_s_at XRCC4 X-ray repair complementing defective repair in 121. 0.2 1.3 0.3 
Chinese hamster cells 4 
208386_x_at DMC1 DMCI1 dosage suppressor of mck1 homolog, —-0.1 1.2 0.4 0.2 
meiosis-specific homologous recombination 
209579_s_at MBD4 Methyl-CpG binding domain protein 4 0.1 0.6 0.2 1 
219317_at POLI Polymerase (DNA directed) iota 0 1 —0.6 0.5 
Other 204178_s_ at RBM14 RNA binding motif protein 14 —0.3 0.5  —0.3 1 
209803_s_at PHLDA2 Pleckstrin homology-like domain, family A, 0.1 0.1 0.4 el 
member 2 
215338_s_at NKTR Natural killer-tumor recognition sequence 0.8 al 0.3 0.2 
Inflammation 
Anti-inflammatory 201502_s_at NFKBIA/IKBA Nuclear factor of kappa light polypeptide gene 0.2 —0.1 1.4 0.5 
enhancer in B-cells inhibitor, alpha/I-kappa B 
alpha 
221577_x_at PLAB Prostate differentiation factor 0.6 —0.2 1.8 0.5 
Proinflammatory 202510_s_at TNFAIP2 Tumor necrosis factor alpha-induced protein 2 -0.3 —0.3 1.2 1.2 
202625_at LYN? v-yes-1 Yamaguchi sarcoma viral related 0.1 —-0.1 —0.2 1.2 
oncogene homology/oncogene LYN 
202637_s at ICAM1? Intercellular adhesion molecule 1 (CD54), human 0.2 —0.3 2.5 1.8 
rhinovirus receptor 
202644_s at TNFAIP3? Tumor necrosis factor alpha-induced protein 3 0.2 —0.1 2.1 1 
203313_s_at TGIF TGFB-induced factor (TALE family homeobox) 0.5 0.3 all pl 
203964_at NMI N-myc (and STAT) interactor (transcription from 0 0.2 0.2 1 
Pol II promoter, inflammatory response, JAK- 
STAT cascade) 
204470_at CXCL1 Chemokine (C-X-C motif) ligand 1 (melanoma 0.8 —0.2 3.1 0.9 
growth-stimulating activity, alpha) 
205114 _s_ at CCL3/MIP1A Chemokine (C-C motif) ligand 3/macrophage 0.5 0.4 14 1 
inflammatory protein 1-alpha 
205476_at CCL20/MIP3A Chemokine (C-C motif) ligand 20/macrophage 0.4 —0.2 1.2 0.8 
inflammatory protein 3-alpha 
205844 at VNNI1 Vanin 1 —0.6 0.3 0.1 12 
205922_at VNN2 Vanin 2 0.7 —0.2 0.8 1 
206336_at CXCL6 Chemokine (C-X-C motif) ligand 6 (granulocyte 0.2 0.4 0.5 13 
chemotactic protein 2) 
207850_at CXCL3/MIP2B Chemokine (C-X-C motif) ligand 3/macrophage 0.2 1.6 5.2 3.5 
inflammatory protein 2-beta 
208991_at STAT3 Signal transducer and activator of transcription 3 0 —0.5 1 0.4 


(acute-phase response factor) 
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Gene 


Description 


Signal log ratio (log, scale)* 
after infection with: 


subcategory Affymetrix HCoV-229E SARS-CoV 
identification Designation 
no. 2h 4h 2h 4h 
209636_at NFKB2” Nuclear factor of kappa light polypeptide gene -0.8 =1.4 2.5 2.7 
enhancer in B cells 2 
209774_x_at CXCL2/MIP2A Chemokine (C-X-C motif) ligand 2/macrophage 1d 0 3 1.6 
inflammatory protein 2-alpha 
209909_s at TGFB2? Transforming growth factor beta 2 0.1 0.2 0.9 21 
211000_s at IL-6ST’ Interleukin-6 signal transducer (gp130, oncostatin 0.5 0.9 0 13 
M receptor) 
211506_s at IL-8? Interleukin-8 0.6 0.7 3.3 21 
212063_at CD44 CD44 antigen (homing function and Indian blood 0.1 —0.1 0.5 11 
group system) 
212989_at CXCL10 Chemokine (C-X-C motif) ligand 10 or Mob 0.2 0.6 0.9 1 
protein 
215101_s_ at CXCLS Chemokine (C-X-C motif) ligand 5 —0.2 0.4 1.6 2 
217738_at PBEF” Pre-B-cell colony-enhancing factor 0 —0.1 AS: 13. 
AFFX- ISGF-3 Interferon-stimulated transcription factor 3 0.1 0.4 0.4 13. 
HUMIS 
GF3A/ 
M97935 
_S_at 
Immuno- 201667_at GJA1/CX43 Gap junction protein, alpha 1/gap junction 0 0.3 0.7 1.4 
modulatory protein 43-KD/connexin 43/heart connexin 
220954 _s at PILR, BETA Paired immunoglobulin-like receptor beta —0.2 1 =03 0.9 
Other 202531_at IRF1 Interferon regulatory factor 1 0.1 —0.6 1 0.1 
204622_x_ at NR4A2 Nuclear receptor subfamily 4, group A, member 2 0:3: =-0.3 =1.6 0.9 
204897_at PTGER4 Prostaglandin E receptor 4 (subtype EP4) 0.2 0.2 12 alu 
205220_at Chemokine receptor Putative chemokine receptor 0.2. —0.1 1.7, 1 
HM74 
206785_s_at KLRC2 Killer cell lectin-like receptor subfamily C, 0.3 2.2 3.8 5.2 
member 2 
220491_at HAMP Hepcidin antimicrobial peptide 0.4 —04 2.1 0.2 
Coagulation 
Anticoagulants 209277_at TFPI2? Tissue factor pathway inhibitor 2 1 0.9 2.1 3 
Procoagulants 202446_s_at PLSCR1/MMTRAI1B Phospholipid scramblase 1 0 0.2 0.6 1.2 
201110_s at THBS1? Thrombospondin 1 0.5 0.2 2.3 1 
201693_s at EGR1? Early growth response 1 0.1 —0.9 1 SL 
202627_s at PATI Plasminogen activator inhibitor 1, serpine 1 -04 —0.6 1.3 18 
Metabolism 
Amino acid 201739_at SGK Serum/glucocorticoid regulated kinase 0.2 —0.2 1 0.4 
206669_at GAD1 Glutamate decarboxylase 1 (brain, 67 kDa) 0.1 2.9 2.8 3 
219795_at SLCOA14 Solute carrier family 6 (neurotransmitter =0.1 0.4 11 2 
transporter), member 14 
222018_at NACA Nascent-polypeptide-associated complex alpha 0.3 1 —0.4 0.4 
polypeptide 
Carbohydrate 202464_s_ at PFKFB3 6-Phosphofructo-2-kinase/fructose-2,6- 0.1 0.3 11 1 
biphosphatase 3 
214850_at SMAS H. sapiens beta glucuronidase pseudogene dt 135 0.8 21 
215977_x_at GK Glycerol kinase 0 0.6 0.8 1 
218918 at MANIC1 Mannosidase, alpha, class 1C, member 1 —0.3 0.2 -08 =1 
218985_at SLC2A8 Solute carrier family 2 (facilitated glucose -0.1 0.4 0.3 1 
transporter), member 8 
219508_at GCNT3 Glucosaminyl (N-acetyl) transferase 3, mucin type 0.2 0.1 1 i 
219797_at MGAT4A Mannosyl (alpha-1,3-)-glycoprotein beta-1,4-N- 0.4 0.7 0.9 1.2 


acetylglucosaminyltransferase, isoenzyme A 
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Gene 


Description 


Signal log ratio (log, scale)* 
after infection with: 


subcategory Affymetrix HCoV-229E SARS-CoV 
identification Designation 
no. 2h 4h 2h 4h 
Lipid 202067_s_ at LDLR? Low-density lipoprotein receptor (familial 0.5 —0.3 15 -0.3 
hypercholesterolemia) 
204242 s at ACOX3 Acyl coenzyme A oxidase 3, pristanoyl 0.3 0 0.5 13. 
205769_at SLC27A2 Solute carrier family 27 (fatty acid transporter), 0 -0.3 -06 =1.2 
member 2 
205822_s_at HMGCS1 3-Hydroxy-3-methylglutaryl coenzyme A synthase 0 0.8 0.7 ll 
1 (soluble) 
212226 _s at PPAP2B? Phosphatidic acid phosphatase type 2B 0 —0.3 12 -0.3 
Transporters 202437_s at CYP1B1°? Cytochrome P450, family 1, subfamily B, 0.6 AS: 3 3.5 
polypeptide 1 
204685_s_at ATP2B2 ATPase, Ca** transporting, plasma membrane 2 0 =03 OL =f 
205195_at AP1S1° Adaptor-related protein complex 1, sigma 1 ber 17 1 21 
subunit 
205749_at CYP1A1 Cytochrome P450, family 1, subfamily A, —0.4 1.7 1.6 2.5 
polypeptide 1 
207604_s at SLC4A7° Solute carrier family 4, sodium bicarbonate 0.3 0.6 0.7 14 
cotransporter, member 7 
208914 at GGA2 Golgi-associated, gamma adaptin ear-containing, 0.2 0.5 0.2 1 
ARF binding protein 2 
210357_s_at SMOX Spermine oxidase 0 —0.2 12 13 
211572_s_ at SLC23A1 Solute carrier family 23 (nucleobase 0.3 0.1 0 bl 
transporters), member 1 
218703_at SEC22L2 SEC22 vesicle trafficking protein-like 2 1 13. =05 0.9 
220786_s_at SLC38A4 Solute carrier family 38, member 4 —0.2 0 -0.7 =1.1 
Miscellaneous 201556_s_at WAMP2 Vesicle-associated membrane protein 2 0 =18 -07 —0:8 
(synaptobrevin 2) 
202238 _s at NNMT° Nicotinamide N-methyltransferase 0.4 0.5 1.9 2.8 
203798_s_at VSNLI1 Visinin-like 1 -0.1 —-O0.1 0 =i 
206293_at SULT2A1 Sulfotransferase family, cytosolic, 2A, —0.2 0 =0:3* sil 
dehydroepiandrosterone preferring, member 1 
208776_at PSMD11 Proteasome (prosome, macropain) 26S subunit, 0 =03. .=02. .=1 
non-ATPase, 11 
210447_at GLUD2 Glutamate dehydrogenase 2 0.2 0.1 09 =1 
213926_s at HRB HIV-1° Rev binding protein 0.2 0.4 0.3 pl 
214837_at ALB Albumin —0.2 1.2 —-0.7 0.2 
216049_at RHOBTB3 Rho-related BTB domain-containing 3 0.5 04  =1 0.6 
Signal transduction 201010_s_at TXNIP Thioredoxin-interacting protein 0.5 0.2 1.6 0.6 
201044 x at DUSP1? Dual-specificity phosphatase 1 0.1 —0.4 3.2 21 
201295_s_at WSB1 SOCS box-containing WD protein SWiP-1 0.2 16 —0.2 14 
201824 at RNF14 Ring finger protein 14 0.2 01  =12 0.3 
203879_at PIK3CD Phosphoinositide-3-kinase, catalytic, delta -0.2 -0.3 0.1 12 
polypeptide 
205205_at RELB v-Rel reticuloendotheliosis viral oncogene 03 =0:5 1.4 0.8 
homolog B, nuclear factor 
of kappa light polypeptide gene enhancer in B 
cells 3 (avian) 
205239_at AREG Amphiregulin (schwannoma-derived growth 0.4 0.8 15 0.7 
factor) 
205868_s at PTPNI11 Protein tyrosine phosphatase, non-receptor type 0.6 1.6 0.6 0.8 
11 (Noonan syndrome 1) 
210056_at RHO6 GTP binding protein 0.7 0.3 3.8 2.6 
211962_s at ZFP36L1 Zinc finger protein 36 0.2 pal 0.7 
212099_at ARHB Ras homolog gene family, member B 0.2 —0.5 1.3 0.6 
212777_at sosi? Son of sevenless homolog 1 (Drosophila) 0 0.8 0.6 12 
213848_at DUSP7 Dual-specificity phosphatase 7 =0:3 O11 —0.2 
213880_at GPR49 G protein-coupled receptor 49 0.2 05  —0.4 17 
218295_s_ at NUP50 Nucleoporin, 50 kDa 0.2 0.8 0.4 1 
219677_at SSB1 SPRY domain-containing SOCS box protein SSB-1 0.4 0.2 14 0.6 
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after infection with: 


subcategory Affymetrix HCoV-229E SARS-CoV 
identification Designation 
no. 2h 4h 2h 4h 
41660_at CELSR1 Cadherin epidermal growth factor LAG seven- 0.7 0.6 0.4 peal 
pass G-type receptor 1 
AFFX- GPR34 G protein-coupled receptor 34 —2.7 3 0.9 3.7 
HUM 
RGE/ 
M10098_ 
5_at 
Transcription 201169 s at BHLHB2? Basic helix-loop-helix domain containing, class B, 2 0 0 13 18 
201473_at JUNB Jun B proto-oncogene 05° =03 2.4 2.2 
201510_at? ELF3* E74-like factor 3 (Ets domain transcription 0.5 0 19 13 
factor, epithelial specific) 
202601_s_ at HTATSF1/TATSF1 HIV TAT specific factor 1 0.7 14 0.9 11 
202672_s_at ATF3 Activating transcription factor 3 0.6 —0.6 2.2 0.1 
202935_s_at SOX9 SRY (sex determining region Y)-box 9 0.4 —0.3 11 -0.3 
(campomelic dysplasia, autosomal sex reversal) 
203751_x_at JUND? Jun D proto-oncogene 0.5 =1.2 14 —-0.4 
203973_s_at CEBPD CCAAT/enhancer binding protein (C/EBP), delta 0.6 —0.6 2.3 15 
204069_at MEIS1 Meis1, myeloid ecotropic viral integration site 1 -1.1 0.3. —0.7 0.1 
homolog (mouse) 
205187_at SMADS5 Mothers against DPP homolog 5 (Drosophila) 0.9 1 a3 1.4 
SMAD 
205398_s_at SMAD3 SMAD, mothers against DPP homolog 3 0.4 0.1 1 i 
(Drosophila) 
205443_at SNAPC1 Small nuclear RNA activating complex, 0.4 0.6 1 0.8 
polypeptide 1 
205471_s_at DACH1 Dachshund homolog (Drosophila) 1 04 -06 —-0O1 =1 
205975_s_at HOXD1 Homeo box D1 0 —0.3 1 0.7 
206036_s_at REL v-Rel reticuloendotheliosis viral oncogene 0.4 0.2 1 0.1 
homolog (avian) 
206175_x_at ZNF222 Zinc finger protein 222 0.2 1 0.1 0.5 
207147_at DLX2 Distal-less homeo box 2 0.7 =1.8 1.5 1:7 
207558_s_at PITX2 Paired-like homeodomain transcription factor 2 —0.2 11 =1 —0.1 
208129_x_at RUNX1 Runt-related transcription factor 1 (acute myeloid 0.5 7 1 0.4 
leukemia 1; am11 oncogene) 
209291 _at ID4 Inhibitor of DNA binding 4, dominant negative 0.2 -0.2 -06 =1 
helix-loop-helix protein 
209357_at” =CITED2? Cbp/p300-interacting transactivator, with Glu/ 0.1 -01 -14 -02 
Asp-rich carboxy-terminal domain, 2 
209505_at NR2F1 Nuclear receptor subfamily 2, group F, member 1 0.1 0:2. =1 0.3 
209757_s_at MYCN v-Myc myelocytomatosis viral related oncogene, 0.2 -0.3 =11 —-0.7 
neuroblastoma 
derived (avian) 
212641_at HIVEP2 HIV-2 1 enhancer-binding protein 2 0.3. —0.1 1.8 0.2 
213668_s at SOX4? SRY (sex-determining region Y) box 4 0.5 0.7 18 1.3: 
213844 at HOXAS Homeo box A5 —0.1 0.1 pal 1 
213906_at MYBL1 v-Myb myeloblastosis viral oncogene homolog —0.5 04  =1 0.2 
(avian)-like 1 
214639_s at HOXA1 Homeo box Al eaAlel: | = 1k2 1.3 1.9 
215073_s_at NR2F2? Nuclear receptor subfamily 2, group F, member 2 0.1 0 =12 =0:5 
215720_s_ at NFYA Nuclear transcription factor Y, alpha 0.4 1.6 0.8 1.4 
218559_s_ at MAFB v-Maf musculoaponeurotic fibrosarcoma —0.4 0 -14 =13 
oncogene homolog B (avian) 
220444 _at ZNF557 Zinc finger protein 557 2.3 2.2 0.9 2.4 
221773_at ELK3 ETS domain protein —0.3 0.1 —0.1 1 
221883_at PKNOX1 PBX/knotted 1 homeobox 1 0.8 0.9 0.2 17 
36711_at MAFF” v-Maf musculoaponeurotic fibrosarcoma (avian) 0.6 =1 19 0 
oncogene family 
40837_at TLE2 Transducin-like enhancer of split 2 [E(sp1) 0.6 =1 0.5 0.5 
homolog, Drosophila} 
44783_s at HEY1 Hairy/enhancer-of-split related with YRPW motif 1 0 0.3 0.6 1.3: 
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after infection with: 


subcategory Affymetrix HCoV-229E SARS-CoV 
identification Designation 
no. 2h 4h 2h 4h 
Cell cycle 201055_s_at HNRPAO Heterogeneous nuclear ribonucleoprotein AO 0.1 0.6 0.5 1 
202150_s_at NEDD9 Neural precursor cell expressed, developmentally 0 —0.8 11 0.8 
downregulated 9 
202768_at FOSB FBJ murine osteosarcoma viral oncogene -03  —0.4 3.8 —0.1 
homolog B 
202769_at CCNG2? Cyclin G2 0.1 0 =02 15 
203120_at TPS3BP2 Tumor protein p53 binding protein, 2 0.3 0.2 plea 0.4 
204137_at TM7SF1 Transmembrane 7 superfamily member 1 17 1.6 15 2.7 
(upregulated in kidney) 
204159_at CDKN2C Cyclin-dependent kinase inhibitor 2C (p18, =0.2 01 =13 =0.7 
inhibits CDK4) 
205027_s_ at MAP3K8 Mitogen-activated protein kinase kinase kinase 8 —-0.1 0.7 1.9 0.7 
205046_at CENPE Centromere protein E 0.1 1,1. -.-=0;1 0.8 
205425_at HIP1 Huntingtin-interacting protein 1 01 = 0.7 0.9 
207530_s_at CDKN2B Cyclin-dependent kinase inhibitor 2B (p15, 0.4 2.3 2.5 3.9 
inhibits CDK4) 
209680_s_at KIFC1 Kinesin family member Cl 0 0.4 —0.6 1 
210136_at MBP Myelin basic protein 0.4 0.5 11 0.9 
212000_at SFRS14 Arginine/serine-rich 14 splicing factor 0.8 1.2 0.4 0.8 
212769_at TLE3 Transducin-like enhancer of split 3 [E(sp1) 0 0.2 122, 0.7 
homolog, Drosophila} 
221841 s at KLF4 Kruppel-like factor 4 0.4 —0.8 4:3 0 
222303_at ETS2 v-Ets avian erythroblastosis virus E26 oncogene =0.2 0.5 1.2 0.8 
homolog 2 
Cell growth and 201329_s at ETS2° v-Ets crythroblastosis virus E26 oncogene 0.3 0.1 2 13 
differentiation homolog 2 (avian) 
201341_at? ENCI/PIG10/NRPB  Ectodermal-neural cortex (with BTB-like domain) 0 0.2 0.8 12. 
p53-induced gene 10; nuclear restricted 
protein/brain 
207401_at PROX1 Prospero-related homeobox 1 0.3 0.5 72 11 
210001_s_ at SOCS1 Suppressor of cytokine signaling 1 0.2 —0.5 1.1 0.4 
210587_at INHBE Inhibin, beta E 0.2 —0.1 0.9 11 
213004_at ANGPTL2 Angiopoietin-like 2 -0.1 0.1 0.3 i 
216684_s at SS18? Synovial sarcoma translocation, chromosome 18 0.9 15 13 17 
Stress response 200664_s_at DNAJB1 Heat shock 40-kDa protein 1 0.1 0.4 1 0.7 
202376_at AACT/SERPINA3 Alpha-1 antichymotrpysin/serine (or cysteine) -04 —-0.8 14 1.3, 
proteinase inhibitor, clade A (alpha-1 
antiproteinase, antitrypsin), member 3 
202581_at HSPA1B/HSP70-2 Heat shock 70-kDa protein 1B/heat shock =0:2 0.8 —0.2 1 
protein, 70 kDa 2 
208607_s_at SAA2 Serum amyloid A2 2.7 2 3.5 3.7 
209799_at PRKAAI Protein kinase, AMP activated, alpha 1 catalytic 0 0.7 0.2 11 
subunit 
215078_at SOD2?’ Superoxide dismutase 2, mitochondrial 0.4 0.4 3 17 
221477_s_at SOD7 Superoxide dismutase 2, mitochondrial 0 0.1 11 1:3, 
Miscellaneous 202241_at TRIB1 Tribbles homolog 1 (Drosophila) -0.1 —0.6 Ad. 0.3 
202458_at SPUVE Protease, serine, 23 0.1 0.2 0.4 1 
202760_s_at AKAP2 A kinase (PRKA) anchor protein 2 0.6 0 pial 0.1 
203234_at UPP1 Uridine phosphorylase 1 —0.6 0.3 0.3 13. 
205890_s_at UBD Ubiquitin D 0 —0.2 a3 2 
206374_at DUSP8 Dual-specificity phosphatase 8 0.6 —0.6 7: ~~ =0:7 
210159 s at TRIM31? Tripartite motif-containing 31 2.2 11 2.3 3 
212027_at RBM25 RNA binding motif protein 25 0.2 1d 0.3 0.9 
212244 at GRINLIA Glutamate receptor, ionotropic, N-methyl D- 0.2 01 =1 0.2 
aspartate-like 1A 
212905_at CSTF2T Cleavage stimulation factor, 3’ pre-RNA, subunit —0.5 0 =i 0.1 


2, 64 kDa, tau variant 
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subcategory Affymetrix HCoV-229E SARS-CoV 
identification Designation 
no. 2h 4h 2h 4h 
213006_at KIAA0146 KIAA0146 protein 0.2 —0.6 1.8 ile 
217475_s_at LIMK2 LIM domain kinase 2 0.2 0.3 1 0.5 
218737_at SBNO1 Sno, strawberry notch homolog 1 0.9 1 0.3 1 
218810_at FLJ23231 Hypothetical protein FLJ23231 = —3.5 2.7 1.3 
220992_s_ at Clorf25 Chromosome 1 open reading frame 25 0.1 12 —-0.1 1 
221258 _s at DKFZP434G2226 Kinesin family member 18A 0 0.6 —0.9 iL 
Unknown 204447_at ProSAPiP1 ProSAPiP1 protein 0.1 —0.8 0 =12 
functions 
204492 _at ARHGAPI1A KIAA0013 gene product 0.5 1 =02 0.5 
208154_at LOCS51336 Mesenchymal stem cell protein DSCD28 0.7 0.3 O14 =1.2 
212686_at PPM1H Protein phosphatase 1H (PP2C domain 0 -—04 -O1  =1.1 
containing) 
212845_at SAMD4 Sterile alpha motif domain containing 4 0.1 —0.4 1 —0.5 
213069_at HEG? HEG homolog 0.3 —-0.3 11 -0.2 
213256_at MGC48332 Hypothetical protein MGC48332 0 0.2 0.4 peal 
214077_x_at MEIS3 Meis1, myeloid ecotropic viral integration site 1 0.5 0.4 0.3 12 


homolog 3 (mouse) 


214862_x_at MRNA; cDNA DKFZp564G1162 0.1 01 =1 0.6 
215281_x_at POGZ Pogo transposable element with ZNF domain 0 1 0.1 —-0.3 
216715_at 1.2 1.3 1 1.1 
217047_s at FAM13A1? Family with sequence similarity 13, member Al —0.3 06 =1 =0;1 
217300_at -1 0.2 =1 0.4 
218031_s at Cl4orf116 Chromosome 14 open reading frame 116 0.1 —0.7 02 =. 
218541_s at CS8orf4 Chromosome 8 open reading frame 4 0.4 0.1 ples) 0.9 
218816_at LANO LAP (leucine-rich repeats and PDZ) and no PDZ 0.1 0.3 0 i 
protein 

219010_at FLJ10901 Hypothetical protein FLJ10901 —0.1 0.1 1 1.6 
219166_at C14orf104 Chromosome 14 open reading frame 104 —0.4 0:3" - =13 0.4 
219696_at FLJ20054 Hypothetical protein FLJ20054 OL =O. =02> 1 
220041_at SMP3 SMP3 mannosyltransferase 04 -11 0.5 =-1 
220441_at FLJ13236 Hypothetical protein FLJ13236 0 0.4 0.5 14 
220770_s_at LOC63920 Transposon-derived Buster3 transposase-like —0.6 13 =0:8 1 
222196_at LOC286434 Hypothetical protein LOC286434 11 0.8 0.9 2.2 
AFFX-r12- Human 18S rRNA? —2.1 AST, 0.7 2.6 

Hs18Sr 

RNA- 

M x at 
AFFX-r2- Human 28S rRNA? 14 0.7 0.5 12 

Hs28Sr 

RNA- 

3_at 


“ Minus, decreased gene expression; no minus, increased gene expression; underlining, two-to fourfold-increased or decreased gene expression; boldface, >4-fold- 


increased or decreased gene expression. 


» More than one probe set has a significant change in gene expression; only the most representative is shown. 


© HIV-1, human immunodeficiency virus type 1. 


tried to study the difference in the expression profile at the 
relatively early stages of the infection of 2 and 4 h postinocu- 
lation, as reported for herpes simplex virus type 1 (17). This 
time frame is biologically relevant because proliferation of the 
Golgi complex and related vesicles and swelling of trans-Golgi 
sacs were observed in infected cells within the first hour of 
infection. Extracellular virus particles were present in 5% and 
30% of the cell populations at 5 and 6 h postinoculation, 
respectively (35). This would also facilitate the analysis, since a 
lower number of altered genes were involved. In the compar- 
ative transcriptomic study, far more genes (n = 136) were 
upregulated by SARS-CoV than by HCoV-229E. In contrast to 
the reported findings of increased antiapoptoptic/inflamma- 


tory gene expression and decreased proapoptotic/inflamma- 
tory gene expression in enterocyte cell lines (10), far more 
proapoptotic and proinflammatory genes were expressed in 
Huh7 cells infected by SARS-CoV but not HCoV-229E. For 
instance, expression of BCL2 was induced by SARS-CoV in 
enterocytes, yet we observed upregulation of its antagonists, 
including BAX and BCL2L11, in Huh7 cells. Moreover, much 
higher expression of other proapoptotic proteins, including 
CASP7, CARD10, PMAIP1, and GADD45B, was also induced 
by SARS-CoV than by HCoV-229E. Furthermore, there was 
marked perturbation of genes involved in cell cycle regulation, 
including induction of the CDKN2B, gene which can mediate 
growth arrest at the G, phase. The induction of proinflamma- 
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FIG. 1. Linkage of the altered genes by using Pathway Assist (version 2.53) shows clustering of genes related to apoptosis, inflammation, and 
coagulation. Only genes with known direct interactions are shown. The software is available from Ariadne Genomics Inc. +, upregulated/increased 


expression; —, downregulated/decreased expression. 


tory cytokines by SARS-CoV was even more prominent when 
compared to HCoV-229E. The induction of IL-8 may be of 
pathogenic importance, as its level has been positively corre- 
lated with disease severity in pulmonary infection by respira- 
tory syncytial virus. Thus, the observed significantly higher 


TFIP2 
PAI1 Coagulation 
THBS1 
IL8 

Inflammation 
NFKB2 
JUNB —— Transcription 
PHLDIA 
CARD10 Apoptosis 
BAX 
PBGD —— House keeping gene 


FIG. 2. Validation of microarray data by semiquantitative RT- 
PCR. The expression patterns of all nine selected genes and the house- 
keeping gene, which were assessed for their relative abundances in 
mock infection, at 2 h and 4 h after SARS-CoV inoculation, and at 2 h 
and 4 h after HCoV-229E inoculation, were consistent with the ex- 
pression pattern shown in the microarray experiment. 


level of IL-8 induced by SARS-CoV compared with HCoV- 
229E in Huh7 cells may recapitulate the host response to these 
viruses by pneumocytes. The induction of various chemokines 
of the CXC or CCL family may mediate the chemotaxis of 
lymphocytes and neutrophils. These alterations in gene expres- 
sion are in keeping with the histological changes of SARS 
hepatitis, in which cellular apoptosis, marked accumulation of 
cells in mitosis with ballooning degeneration of hepatocytes, 
and moderate lymphocytic infiltration were found in liver tis- 
sue biopsy samples (4). Recently, a vaccine study in ferrets 
immunized with the modified vaccinia virus Ankara carrying 
the spike protein showed that the ferrets developed hepatitis 
after being challenged with wild SARS-CoV. This finding tends 
to suggest that the Huh7 cell line model may have some rele- 
vance to the pathogenesis of SARS (50), but our present find- 
ings still may not be directly applicable to the pneumocytes 
involved in SARS. 

The upregulation of genes involved in procoagulation and 
platelet activation was interesting. TFPI2 inhibits thrombin 
generation by binding and inactivation of the tissue factor- 
factor VIIa complex. Upregulation of the gene probably rep- 
resents an inhibitory response to restrain the activation of the 
coagulation pathway during acute inflammation. In contrast, 
TFPI2 also inhibits both free and matrix/cell-associated plas- 
min, thus favoring fibrin deposition, and may have positive role 
in matrix turnover (11). Upregulation of the PAI1 gene ac- 
companied by a dramatic increase in the protein level results in 
an antifibrinolytic response, favoring fibrin deposition during 
the acute inflammatory phase of the disease. It is important to 
remember that mouse hepatitis virus can activate the immune 
coagulation system by the fgl2 gene, encoding a prothrombi- 
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FIG. 3. Fold changes of expression of nine selected genes by RT-qPCR. The P values of the differences in fold changes of the target gene 
expressions are as follows: TFPI2 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.02; TFPI2 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.02; 
PAI (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.05; PAI (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.05; THBS1 (SARS-CoV at 2h 
versus HCoV-229E at 2 h), 0.03; THBS1 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.35; IL-8 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 
<0.01; IL-8 (SARS-CoV at 4 h versus HCoV-229E at 4 h), <0.01; NFKB2 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.05; NFKB2 
(SARS-CoV at 4 h versus HCoV-229E at 4h), <0.01; JUNB (SARS-CoV at 2 h versus HCoV-229E at 2 h), <0.01; JUNB (SARS-CoV at 4 h versus 
HCoV-229E at 4 h), 0.28; PHLDIA (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.20; PHLDIA (SARS-CoV at 4 h versus HCoV-229E at 4 h), 
<0.01; CARD10 (SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.03; CARD10 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.04; BAX 
(SARS-CoV at 2 h versus HCoV-229E at 2 h), 0.04; BAX (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.24. 
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FIG. 4. Validation of protein expression for TFP1, PAI, and IL-8 by enzyme immunoassay. The P values of the differences in protein expression are 
as follows: TFPI2 (SARS-CoV at 2 h versus HCoV-229E at 2 h), <0.01; TFPI2 (SARS-CoV at 4 h versus HCoV-229E at 4 h), 0.37; TFPI2 (SARS-CoV 
at 12 h versus HCoV-229E at 12 h), 0.18; TFPI2 (SARS-CoV at 24 h versus HCoV-229E at 24 h), 0.12; PAI1 (SARS-CoV at 2 h versus HCoV-229E 
at 2 h), <0.01; PAIL (SARS-CoV at 4 h versus HCoV-229E at 4 h), <0.01; PAIL (SARS-CoV at 12 h versus HCoV-229E at 12 h), <0.01; PAIL 
(SARS-CoV 24 at h versus HCoV-229E at 24 h), 0.01; IL-8 (SARS-CoV at 2 h versus HCoV-229E at 2 h), <0.01; IL-8 (SARS-CoV at 4 h versus 
HCoV-229E at 4 h), <0.01; IL-8 (SARS-CoV at 12 h versus HCoV-229E at 12 h), <0.01; IL-8 (SARS-CoV at 24 h versus HCoV-229E at 24 h), <0.01). 
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nase (12). This enzyme can induce macrophage procoagulation 
activity which results in fibrin deposition on the endothelia of 
intrahepatic veins and hepatic sinusoids. The end result could 
be confluent hepatocellular necrosis. The low number of liver 
biopsies done with these patients may account for the lack of 
reports on these full-blown changes related to vascular dam- 
age. However, systemic vasculitis, including edema, localized 
fibrinoid necrosis, and infiltration by monocytes, lymphocytes, 
and plasma cells into vessel walls in various tissues, was re- 
ported (13). Thrombosis was found in small veins. Interest- 
ingly, marked upregulation of a proapoptotic gene, PHLDA1, 
was observed in SARS-CoV infection of Huh7. A previous 
study has shown that overexpression of this gene in vascular 
endothelial cells would lead to decreased cell adhesion and 
induce detachment-mediated apoptosis. This gene, if similarly 
induced in vascular endothelial cells infected by SARS-CoV, 
may contribute to the vascular damage induced by SARS-CoV 
infection (21). 

The clinical manifestations and histological changes involv- 
ing the pneumocytes, enterocytes, and vascular endothelium 
were not surprising because ACE2, the host cell receptor for 
viral entry, was present in all these cell types (19). Similarly, 
abnormal urinalysis associated with increased viral load in 
urine is expected, since SARS-CoV could be cultured from or 
antigen expression could be detected in kidney cell lines and 
kidney tissues of patients (22). However, neither ACE2 nor 
L-SIGN could be detected by immunohistochemical staining in 
normal human hepatocytes or Huh7 (1, 42). This could be 
explained by an alternative cellular receptor for viral entry or 
just because the amount of receptor required for cell entry is so 
low that the cells may appear negative on immunohistochem- 
ical staining. Our previous study has shown that a high viral 
load in serum or stool is associated with hepatic dysfunction 
(22). The viral load in serum reflects the dynamic of the viral 
replication from any organs and the clearance mechanism of 
the host, whereas a high viral load in the stool may indirectly 
reflect the degree of portal venous viremia. Thus, the lack of 
cell death and inflammation in the intestine may still have 
importance in the pathogenesis of SARS in other body sites. 

In conclusion, SARS-CoV produces more severe distur- 
bance of host cell gene expression in a human epithelial cell 
line of liver origin than does HCoV-229E at the early stage of 
infection. There are marked alterations in gene expression 
related to apoptosis, inflammation, and procoagulation. These 
findings are consistent with the histological changes of SARS, 
especially in the liver and blood vessels. Besides antivirals (6, 
24), other modalities of treatment, such as antiapoptotic agents 
(34), immunomodulators against inflammation (7), and modi- 
fiers of coagulation, should be considered in future research. It 
is important to remember that many patients continued to 
deteriorate despite a decreasing viral load 2 to 3 weeks after 
the onset of SARS. 
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